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This article describes the search for pair production of scalar leptoquarks performed by the CMS 
collaboration using the data from the 2010 proton-proton run of the Large Hadron Collider at a 
center of mass energy of 7 TeV. The following final state signatures from the decay of a pair of 
scalar leptoquarks are described in here: the di-lepton plus jets final state (where the leptons are 
either both electrons or both muons) and the lepton plus jets and transverse missing energy final 
state (where the lepton is an electron). 

I. INTRODUCTION 

Leptoquarks (LQ) are hypothetical particles that are predicted by many extensions [TJJ5] of the Standard 
Model of particle physics (SM), such as Grand Unification Theory, Technicolor and Composite models. They 
carry both baryon and lepton numbers and thus couple to a lepton and a quark. They carry fractional electric 
charge, they are color triples, and can have either zero or one unit of spin (i.e. can be either scalar or vector 
particles). Existing experimental limits on lepton number violation, Flavor Changing Neutral Current, proton 
decay, and other rare processes favor three generations of leptoquarks, with no inter-generational mixing. 

The production and decay of LQs are characterized by: the mass of the LQ particle (M^q), its decay branching 
ratio into a lepton and a quark (usually denoted as /3) and the Yukawa coupling at the LQ-lepton-quark vertex 
(A). At hadron colliders, leptoquarks are mainly produced in pairs, via gluon-gluon fusion and quark anti-quark 
annihilation. The dominant pair production mechanisms do not depend strongly on A and the single production 
of leptoquarks does not become significant (and it does not in any case invalidate the search results presented 
here) in the range of LQ masses probed with the 2010 LHC data. 

The final state event signatures from the decay of pair produced LQs can be classified as: dilepton and jets 
(both LQ and anti-LQ decay into a charged lepton and a quark); lepton, missing transverse energy and a jet (one 
LQ decays into a charged lepton and a quark, while the other decays into a neutrino and a quark) ; and missing 
transverse energy and jets (both LQ and anti-LQ decay into neutrinos and quarks). The three signatures have 
branching ratios of /3 2 , 2/3(1-/3), and (1 — /3) 2 , respectively. The charged leptons can be either electrons, muons, 
or tau leptons, corresponding to the three generation of LQs. Only electrons and muons are considered here. 

II. SEARCH FOR PAIR-PRODUCTION OF SCALAR LEPTOQUARKS WITH THE CMS 

DETECTOR 

Searches for first and second generation pair-produced scalar LQs were performed and published on 34-36 
pb _1 of 7 TeV proton-proton collision data recorded by the CMS detector [5] during the 2010 LHC run. First 
generation results include searches both in the dilepton and jets final state (eejj), and in the lepton, missing 
transverse energy and jets final state (evjj). Second generation result consist of a search in the dilepton and 
jets final state (/J-fJ-jj)- The eejj and evjj results were combined to attain the best possible exclusion reach in 
all of the parameter space of the first generation pair produced scalar LQ search, i.e. (/3 and Mlq). 

The analysis in all of the channels aims at identifying the existence of new heavy particles by establishing 
an excess of events characteristic of the decay of heavy objects. The analysis starts with using either a single 
or a double lepton trigger path, which is robust and very efficient. An event signature preselection isolates 
events with high transverse momentum final state objects (two or more isolated leptons and two or more jets, 
or one isolated lepton, two or more jets, and significant missing transverse energy indicative of the presence 
of a neutrino). Event variables are identified to effectively separate a possible LQ signal from standard model 
backgrounds (these are the St, M(ll), and Mt{Iv) variables described below) and lower thresholds are placed on 
these variables at preselection level. Backgrounds from standard model sources are estimated and first compared 
with data at the preselection level. Major sources of background are Z/~f*+]ets and TF+jets processes and tt 
(with single top production, diboson processes, and QCD multijet processes being smaller contributions). Major 
backgrounds are either directly determined from data control samples or determined with Monte Carlo samples 
(therefore using kinematic shape information from the MC) but normalized to data in selected control regions. 

A cut-based variable approach is used for the final selection, where the selection is optimized for different LQ 
mass hypotheses by minimizing the expected upper limit on the LQ cross section in the absence of an observed 
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signal using a Bayesian approach [5J [5]. The three variables used in the optimization are: the scalar sum of 
the final objects transverse momenta, St', the invariant mass of the dilepton pair, Mu, for the eejj and /J-fJ-jj 
channels; and the transverse mass, MtQv), of the lepton and neutrino in the evjj analysis. 

After final selection, the data are well described by the standard model background predictions. In the 
absence of an observed signal, an upper limit on the LQ cross section is set using a Bayesian method [H Hj with 
a flat signal prior. A lognormal prior is used to integrate over the nuisance parameters. Using Poisson statistics, 
a 95% confidence level (C.L.) upper limit is obtained on the LQ pair-production cross-section times branching 
ratio as a function of LQ mass. Comparing with the NLO predictions |7] for the scalar LQ pair production 
cross section a lower limit on Mtq is determined for /3 = 1 and f3 = 0.5 (electron only, for the latter). In the 
electron channel, the eejj and evjj channels are combined to further maximize the exclusion in (3 and Mlq, 
specially for the case of /3 = 0.5, where the eejj adds to the sensitivity of the evjj channel. 

III. eejj AND fi/jtjj CHANNELS 

The eejj analysis requires the presence of a single or double electromagnetic trigger (with an efficiency close 
to 100%), two or more isolated electrons with pt > 30 GeV and \t)\ < 2.5, and two or more jets with pr > 30 
GeV and \rj\ < 3.0. In addition, at pre-selection level, Ai?(e,j) > 0.7 is required together with a minimum 
threshold of M ee > 50 GeV and St = Pt(^i) + Pt{^2) + PT(j e h) + PT(jet 2 ) > 250 GeV. Data and standard 
model background predictions agree well at the level of pre-selection, as shown in Fig. [T] and [2j This channel 
uses an integrated luminosity of 33.2 pb^ 1 . 




FIG. 1: Leading electron pr at preselection level for the eejj analysis. 




FIG. 2: St at preselection level for the eejj analysis. 

Similarly, the p/jijj analysis starts from the requirement of a single muon trigger (with an efficiency of 99%), 
two or more isolated muons with pt > 30 GeV and < 2.4 (one of which must be within |f?| < 2.1), and 
two or more jets with pt > 30 GeV and |f?| < 3.0. The two muons are required to be separated in R by at 
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least 0.3. Minimum thresholds of > 50 GeV and St — Pt{hi) + Pt{^2) + Prijetx) + Prijefa) > 250 

GeV complete the preselection requirements. Good agreement is observed between data and standard model 
background predictions at this level (Fig.s^andffl). This channel uses an integrated luminosity of 34.0 pb _1 . 
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FIG. 3: Leading muon pr at preselection level for the nfijj analysis. 



> 




100 200 300 400 500 600 700 800 900 1000 



S T (GeV) 

FIG. 4: St at preselection level for the fJ.fJ.jj analysis. 

The major backgrounds from standard model processes to the dilepton and jets channels come from Z/j* + 
jets and it production. The Z/j* + contribution is determined from MC, but it is rescaled to the data in the 
Z peak region at preselection. The invariant mass of the dilepton pair is shown in Fig. [5] and in Fig. [6] for 
the eejj and fJfjjj channel, respectively. The normalization factors are 1.20 ± 0.14 for eejj and 1.28 ± 0.14 
for fJ.fJ.jj- The systematic uncertainty on the Z/j*+ jets background prediction in both channels comes from 
the statistically dominated uncertainty in the normalization factors and from a shape uncertainty obtained 
by comparing the yields of Z/~/*+ jets MC samples generated with different renormalization and factorization 
scales and matching thresholds. The it normalization uncertainty is based on a CMS measurement of the it 
cross-section [10]. Smaller W+ jets and diboson plus jets backgrounds are determined entirely from MC and are 
found to be negligible. The background contributions from QCD multijet processes are determined from data 
control regions and are also found to be negligible. At final selection, a dilepton invariant mass cut (M ee >125 
GeV and M w >115 GeV) is placed to suppress most of the Z/-f* + jets background and a St threshold is 
optimized for each M^q hypothesis. St is effective at removing most of the it background and any residual 
backgrounds surviving the preselection criteria and the cut on the dilepton invariant mass. The number of data, 
signal MC, and background events for each of the final selection criteria optimized for different Mlq are listed 
in Table [JJ and Table [TT] for the eejj and fJfijj channels respectively, together with the optimized St thresholds. 

In absence of an excess above standard model backgrounds expectation, an upper limit on the LQ cross section 
is set using a Bayesian method [HI [H] with a flat signal prior. A log-normal probability density function is used 
to integrate over the systematic uncertainties. Major sources of systematic uncertainties for both channels are 
the uncertainties on the determination of the luminosity, the jet energy scale, the dilepton reconstruction and 
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FIG. 5: Invariant mass of the two highest pr electrons at preselection level for the eejj analysis. 




FIG. 6: Invariant mass of the two highest pr muons at preselection level for the analysis. 



identification efficiencies, and the uncertainties associated with the normalization and modeling of the main 
backgrounds, Z/7*+jets and it. Using Poisson statistics, a 95% confidence level (C.L.) upper limit is obtained 
on a x f3 2 . This is shown in Fig. [7] and Fig. [8] together with the NLO [7] predictions for the scalar LQ pair 
production cross section. The systematic uncertainties are included in the calculation as nuisance parameters. 
With the assumption that /3 = 1, first generation and second-generation scalar leptoquarks with masses less 
than 384 and 394 GeV are excluded at 95% C.L. [UlIU]. This is in agreement with the expected limits of 391 
and 394 GeV. 



A. evjj channel and combination with eejj 



The evjj analysis requires the presence of a single electromagnetic trigger, one isolated electrons with pr > 35 
GeV and \ij\ < 2.2, and two or more jets with px > 30 GeV and \r}\ < 3.0. In addition, at pre-selection level, 
AR(e,j) > 0.7 is required together with transverse missing energy, MET > 50 GeV, a veto on the presence of 
muons in the event, St = Pr{e) + M ET +pr(jeti) +pr(jet2) > 250 GeV, and azimuthal cuts between final state 
objects, \A<p(MET, e)| > 0.8 and \A<p(MET,jeti)\ > 0.5 to reduce the contribution from mis-reconstructed 
events. Data and standard model background predictions agree well at the level of pre-selection, as shown in 
Fig. [9] This channel uses an integrated luminosity of 36 pb _1 . 

The major backgrounds from standard model processes to the evjj channel are W+ jets and ti production. 
The VF+jets contribution is determined from MC and rescaled to the data in the 50< Mt(cv) <110 GeV 



region at preselection. The Mt{cv) distribution is shown in Fig. 10 The systematic uncertainty on the W+ 
jets background prediction comes from the statistically dominated uncertainty in the normalization factors and 
from a shape uncertainty obtained by comparing the yields of W+ jets MC samples generated with different 
renormalization and factorization scales and matching thresholds. The ti normalization uncertainty is based 
on a CMS measurement of the tt cross-section [10] and the shape uncertainty is determined from MC. Other 
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TABLE I: Number of eejj events for LQ signal, backgrounds, and data after full selection with 33.2 pb -1 . The product 
of signal acceptance and efficiency is also reported for different LQ masses. The Z/7*+jets MC has been normalized to 
the data as described in the text. Other backgrounds include W + jets, di-boson, and single top. The uncertainties are 
statistical. The observed and expected 95% C.L. upper limit (u.l.) on the LQ pair production cross section a are shown 
in the last column. 







Signal Samples (MC) 


Standard Model Background Samples (MC) 


Events 


Obs./Exp. 


(S T Cut) 


Selected 


Acceptance 




Selected Events in 




in 


95% C.L. 




[GeV] 


Events 


X Efficiency 


tt + jets 


Zh"+ jots 


Others 


All 


Data 


u.l. on a [pb] 


200 


(S T >340) 


117.5±0.8 


0.297±0.002 


2.6 ±0.1 


2.0 ±0.2 


0.27±0.05 


4.9 ±0.2 


2 


0.441 / 0.720 


250 


(S t >400) 


43.8±0.2 


0.380±0.002 


1.3 ±0.1 


1.3 ±0.1 


0.14±0.02 


2.7 ±0.1 


1 


0.309 / 0.454 


280 


(S T >450) 


24.4±0.1 


0.403±0.002 


0.69±0.05 


0.87±0.07 


0.10±0.02 


1.7 ±0.1 


1 


0.305 / 0.373 


300 


(S T >470) 


17.3±0.09 


0.430±0.002 


0.52±0.05 


0.75±0.07 


0.10±0.02 


1.4 ±0.1 


1 


0.292 / 0.332 


320 


(St >490) 


12.3±0.06 


0.451±0.002 


0.43±0.04 


0.65±0.07 


0.08±0.02 


1.2 ±0.1 


1 


0.283 / 0.305 


340 


(St >510) 


8.88±0.04 


0.469±0.002 


0.32±0.04 


0.56±0.06 


0.08±0.02 


0.96±0.08 


1 


0.278 / 0.279 


370 


(St >540) 


5.55±0.02 


0.496±0.002 


0.26±0.03 


0.47±0.06 


0.07±0.02 


0.80±0.07 


1 


0.267 / 0.254 


400 


(S t >560) 


3.55±0.02 


0.522±0.002 


0.20±0.03 


0.41±0.05 


0.06±0.02 


0.67±0.07 


1 


0.257 / 0.234 


450 


(St >620) 


1.70±0.01 


0.539±0.002 


0.12±0.02 


0.28±0.05 


0.02±0.01 


0.42±0.06 





0.174 / 0.210 


500 


(St >660) 


0.868±0.003 


0.565±0.002 


0.08±0.02 


0.23±0.05 


0.02±0.01 


0.33±0.05 





0.166 / 0.194 



TABLE II: The fifijj data event yields in 34.0 pb" 1 for different LQ mass hypotheses, together with the optimized St 
threshold (in GeV) for each mass, background predictions, expected LQ signal events (S), and signal selection efficiency 
times acceptance (es). Mlq and St values are listed in GeV. The Z/j* — > nn + jets and tt contributions are rescaled by 
the normalization factors determined from data. Other backgrounds correspond to VV, W +jets, and multijet processes. 
Uncertainties are statistical. 





M LQ 


MC Signal Samples 


Monte Carlo Background Samples 


Events 


Obs./Exp. 


(S T Cut) 


Selected 


Acceptance 




Selected Events in 




in 


95% C.L. 




[GeV] 


Events 


X Efficiency 


tt -f- jets 


Z/ 7 * +jets 


Others 


All 


Data 


u.l. on a [pb] 


200 


(St > 310) 


160±20 


0.388±0.003 


4.6±0.1 


4.08±0.07 


0.1±0.01 


8.8±0.2 


5 


0.438 / 0.695 


225 


(S T > 350) 


89±9 


0.421±0.003 


3.1±0.1 


2.99±0.05 


0.07±0.01 


6.2±0.1 


3 


0.339 / 0.547 


250 


(S T > 400) 


51±5 


0.437±0.003 


1.88±0.09 


1.92±0.04 


0.051±0.009 


3.9±0.1 


3 


0.366 / 0.436 


280 


(S T > 440) 


28±3 


0.467±0.003 


1.15±0.07 


1.53±0.03 


0.038±0.008 


2.72±0.08 


3 


0.371 / 0.361 


300 


(S T > 440) 


21±2 


0.518±0.004 


1.15±0.07 


1.53±0.03 


0.038±0.008 


2.72±0.08 


3 


0.335 / 0.326 


320 


(S T > 490) 


14±1 


0.509±0.004 


0.64±0.05 


1.12±0.02 


0.019±0.005 


1.78±0.06 


2 


0.300 / 0.292 


340 


(S T > 530) 


9±1 


0.508±0.003 


0.4±0.04 


0.79±0.01 


0.01±0.004 


1.20±0.04 


1 


0.245 /0.264 


400 


(St > 560) 


4.0±0.4 


0.578±0.004 


0.31±0.04 


0.67±0.01 


0.01±0.004 


0.99±0.04 


1 


0.219 / 0.222 


450 


(St > 620) 


1.9±0.2 


0.600±0.004 


0.19±0.03 


0.49±0.01 


0.006±0.003 


0.69±0.03 





0.153 / 0.199 


500 


(St > 700) 


0.9±0.1 


0.602±0.004 


0.09±0.02 


0.277±0.006 


0.003±0.002 


0.37±0.02 





0.152 / 0.180 



sources of background determined from MC are found to be negligible. The background contributions from 
QCD multijet processes are determined from data control regions and are also found to be negligible. 

At final selection, the transverse mass of the electron and neutrino Mt{cv) is required to be above 125 GeV 
and a min(j>r(e), MET) >85 GeV cut is placed to suppress most of the W+ jets background. A St threshold 
is then optimized for each Mlq hypothesis. St is effective at removing most of the tt background and any 
residual backgrounds surviving the preselection criteria and the W veto cuts listed above. The number of data, 
signal MC, and background events for each of the final selection criteria optimized for different Mtq are listed 
in Table [TTTJ together with the optimized St thresholds. 

In absence of an excess above standard model backgrounds expectation, an upper limit on the LQ cross 
section is set using a Bayesian method similar to the one used in the eejj analysis. The evjj channel has 
maximum sensitivity for /3 = 0.5, and also offers sensitivity for lower values of (3, where the eejj channel 
quickly runs out of sensitivity due to the quadratic behavior of the branching fraction of the LQ pair into the 
above channel. Thus, combining it with the eejj channel allows one to extend the sensitivity of the search in 
the intermediate /3 range compared to either of the two analysis considered separately. The two channels are 
combined using the same Bayesian approach used to set the individual limits. Assuming a flat prior for the 
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FIG. 7: The expected and observed 95% C.L. upper limit on the first-generation scalar LQ pair production a x ft 2 as 
a function of the LQ mass together with the NLO theoretical cross section curve. The shaded band on the theoretical 
values includes CTEQ6 PDF uncertainties and the error on the LQ production cross section due to re-normalization and 
refactorization scale variation. The results correspond to 33.2 pb~ in the eejj channel. 
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FIG. 8: The expected and observed 95% C.L. upper limit on the second-generation scalar LQ pair production a x f3 2 as 
a function of the LQ mass together with the NLO theoretical cross section curve. The shaded band on the theoretical 
values includes CTEQ6 PDF uncertainties and the error on the LQ production cross section due to re-normalization and 
refactorization scale variation. The results correspond to 34.0 pb" 1 in the W-ijj channel. 



signal cross section, the expected signal yield in the eejj [evjj] channel is found by multiplying the cross section 
by the signal efficiency, integrated luminosity, and the branching fraction /3 2 [2/3(1-/3)]. Since the major sources 
of systematic uncertainty between the two channels arc correlated, a combined likelihood is constructed from 
the individual likelihoods for the two channels. The 95% C.L. mass limit as a function of /3 is shown in Fig. [lTj 
along with the individual limits from the two channels. With the assumption that /3 = 0.5, first generation 
scalar leptoquarks with masses less than 340 GeV are excluded at 95% C.L. [T3"] . 
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FIG. 9: Missing transverse energy at preselection level for the ei/jj analysis. 




m T ,ev [GeV] 



FIG. 10: Transverse mass of the electron and neutrino, Mr(eu), at preselection level for the eujj analysis. 



B. Conclusions and outlook 



Searches for pair production of first and second generation scalar leptoquarks have been performed at CMS 
in the final states with: two charged leptons and two jets, one charged lepton (electron), missing transverse 
energy and two jets, with the full 2010 statistics. In the absence of a signal we exclude first generation LQ with 



TABLE III: Number of eujj events for the first generation LQ signal, backgrounds, and data samples after the full 
analysis selection. The optimum St threshold is reported for each LQ mass. All uncertainties are statistical. The 
product of signal acceptance and efficiency is also reported for different LQ masses. 



A/lq 


MC Signal Samples 


MC and QCD Background Samples 


Events 


(S t cut) 


Selected 


Acceptance 




Selected Events in 




in 


[GeV] 


Events 


X Efficiency 


tt -+- jets 


iy+jcts 


Other Bkgs 


QCD 


All Bkgs 


Data 


200 (S T > 350) 


34.5±0.2 


0.161 


3.6±0.1 


2.2±0.3 


0.48±0.06 


0.20±0.04 


6.5±0.3 


5 


250 (S T > 410) 


15.9±0.1 


0.255 


2.24±0.09 


1.7±0.3 


0.35±0.05 


0.18±0.05 


4.4±0.3 


3 


280 (S T > 460) 


9.54±0.05 


0.291 


1.43±0.08 


1.2±0.2 


0.29±0.05 


0.14±0.04 


3.1±0.2 


3 


300 (S T > 490) 


6.89±0.03 


0.317 


1.09±0.07 


1.0±0.2 


0.27±0.05 


0.14±0.04 


2.5±0.2 


2 


320 (S T > 520) 


5.03±0.02 


0.339 


0.75±0.05 


0.8±0.2 


0.22±0.05 


0.13±0.04 


1.9±0.2 


2 


340 (S T > 540) 


3.73±0.02 


0.364 


0.65±0.05 


0.7±0.2 


0.20±0.05 


0.12±0.04 


1.6±0.2 


2 


370 (S T > 570) 


2.40±0.01 


0.396 


0.50±0.04 


0.6±0.1 


0.18±0.04 


0.08±0.03 


1.3±0.2 


1 


400 (S T > 600) 


1.57±0.01 


0.426 


0.34±0.04 


0.5±0.1 


0.17±0.04 


0.08±0.03 


l.liO.l 


1 


450 (S T > 640) 


0.797±0.003 


0.467 


0.26±0.03 


0.4±0.1 


0.13±0.04 


0.08±0.04 


0.9±0.1 





500 (S T > 670) 


0.417±0.001 


0.500 


0.18±0.03 


0.4±0.1 


0.12±0.04 


0.08±0.04 


0.8±0.1 
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FIG. 11: The combined expected and observed 95% C.L. mass limit for first-generation scalar leptoquarks as a function 
of p. Individual channel limits (eejj and evjj) are also shown. 

masses below 384 GeV (0 = 1), second generation LQ with masses below 394 GeV (/? = 1), and first generation 
LQ with masses below 340 GeV ((3 = 0.5). The j3 = 1 results were the most stringent limits at the time of 
publication. The f3 = 0.5 and f3 = 1 results have been combined and have been submitted for publication. The 
analysis of the 2011 data is ongoing. 



[1] S. Dimopoulos et al, "Mass Without Scalars", Nucl. Phys., B155 (1979). 
[2] S. Dimopoulos, " Technicolored signatures", Nucl. Phys, B168 (1980). 

[3] E. Eichten et al, "Dynamical breaking of the weak interaction symmetries", Phys. Lett., B90 (1980). 
[4] V.D. Angelopoulos el al, "Search for new quarks suggested by the superstring" , Nucl. Phys., B292 (1987). 
[5] W. Buchmiiller, W. el al, "Leptoquarks in Lepton-Quark Collisions", Phys. Lett., B191 (1987). 
[6] R. Adolphi et al, "The CMS experiment at the CERN LHC", JINST, 0803 (2008). 
[7] M. Kramer et al, "Pair production of scalar leptoquarks at the CERN LHC", Phys. Rev. D, 71 (2005). 
[8] Particle Data Group, "Review of particle physics - Section 32.3.1: The Bayesian approach", Phys. Lett., B667 
(2008). 

[9] I. Bertram et al, "A Recipe for the Construction of Confidence Limits", FERMILAB-TM-2104 (2000). 
[10] V. Khachatryan et al, CMS Collaboration, "First measurement of the cross section for top-quark pair production in 

proton-proton collisions at yfs = 7 TeV", Phys. Lett. B, 695 (2011). 
[11] V. Khachatryan et al, CMS Collaboration, "Search for pair production of first-generation scalar leptoquarks in pp 

collisions at s/s = 7 TeV", Phys. Rev. Lett., 106 (2011). 
[12] V. Khachatryan et al, CMS Collaboration, "Search for pair production of second-generation scalar leptoquarks in 

pp collisions at y/s = 7 TeV", Phys. Rev. Lett., 106 (2011). 
[13] S. Chatrchyan et al, CMS Collaboration, "Search for first generation scalar leptoquarks in the evjj channel in pp 

collisions at yfs = 7 TeV", Phys. Rev. B, 703 (2011). 




